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We investigate the stability of boron fuUerene sets Brg, Byg and B82- We evaluate the ground 
state energies, nucleus-independent chemical shift (NICS), the binding energies per atom and the 
band gap values by means of first-principles methods. We construct our fuUerene design by capping 
of pentagons and hexagons of Beo cage in such a way that the total number of atoms is preserved. In 
doing so, a new hole density definition is proposed such that each member of a fuUerene group has 
a different hole density which depends on the capping process. Our analysis reveal that each boron 
fuUerene set has its lowest-energy configuration around the same normalized hole density and the 
most stable cages are found in the fuUerene groups which have relatively large difference between 
the maximum and the minimum hole densities. The result is a new stability measure relating the 
cage geometry characterized by the hole density to the relative energy. 

The geometric structure of boron allotropes in respect to their ground state energy is a subject of great interest. 
Studies focused on searching the lowest energy configurations has yielded a variety of stable boron structures due to 
electron deficient nature of boron. Such structures involve molecular wheels [DH], boron sheets composed of triangular 
and hexagonal motifs [SHU] , cage architectures based on empty polygons [13] , crossing double rings [HI [H] , snub 
structures [8|, filled hexagons [1^ and filled pentagons [171 EH]- Boron based forms are known for their complex 
chemistry. Unlike carbon, both pure boron and many boron rich compounds contain B12 cage structures in their 
crystalline state. Furthermore, boron crystals show peculiarities such as uncommon 2c-le and 3c-2e bonding [IS] 
and unique aromaticity [IHIIH]- Most crystal samples exhibit defects, vacancies and interstitials, due to the flexible 
bond structure [22]. Among such structures, the newest form, orthorhombic 7-B28, is found to be the second hardest 
elemental material after diamond with a value of 58 GPa [231 HI] and has strong covalent interatomic interactions 
even higher than the intraicosahedral bonds |2S] . 

A key concept in stability analysis of boron sheets is the hole density. This is defined as the ratio of the number 
of hexagonal holes to the number of atoms in the triangular sheet. The so called a-sheet was found to be the most 
stable boron sheet with 77= 1/9 ^TU]. In a recent work, it is shown that hexagonal holes of a-sheet are the scavengers 
of extra electrons from the filled hexagons [26) : a unique aspect of a-sheet, in conjunction with the hole density, is 
that the ratio 77=1/9 exactly corresponds to the ratio of the extra tt electron to the number of a electrons in a filled 
hexagon. Boron layers composed of hexagonal hole doped triangular lattices manifest polymorphism |27j i.e. multiple 
2D boron layers have comparable stabilities even some of them have lower energy than the a sheet [251 US] • Boron 




B76'°^'(C3) Brs'^-'{C^) B82"^'(C2;,) 

5.73 5.76 5.77 



FIG. 1: (Color online) Optimized structures with Et values using PBE/6-31(d; p) method of energetically most stable isomers 
a) B^r(C.), b) B^i^«(Ci) and c) Bj|^«(C2h), respectively 
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sheets g(l/8), g(2/15), ai and 012 are the examples of such ground state structures. Such correlation between the 
hole density and energy is a generic feature of 2D boron sheets. Different types of correlations may exist in different 
boron structures. 3D crystal forms have correlations between band structure and defects [55]. In small boron clusters, 
aromacity and planarity are strongly related j30]. In this respect, we aim to find such correlations in boron cages. 

Stability analysis of boron fuUerenes, on the other hand, is a formidable task. Szwacki et al. showed that B8o(I/i) 
buckyball composed of boron is stable [15] . Its structure contains 20 additional capped atoms on the hexagons of Bgo 
which is a structural analog of Cgo- Different structural patterns of Bgo might give lower energy cages like volleyball 
shaped Bgo [17) and other Bgo isomers i5T| . Recently, Pochet et al. has found an energy comparability between 
ordered and disordered arrangements in Bgg |32j which is, in principle, an analog of the polymorphism encountered 
in the boron sheets. The structural properties of fullerenes are intimately related to those of sheets as the latter can 
be viewed as the building blocks of the former. Nonetheless, not all the key properties of boron fullerene structures 
may be derived from boron sheet structures for the most stable boron sheets could not form the most stable boron 
cages [H]. In this respect, stability analysis of boron fullerenes requires a proper definition of hole density for the 
cage structures which will be introduced shortly. 

In our analysis we present three groups of boron fullerenes, By^-^, B^-^ and B^-^, where X and Y values are the 
number of filled hexagons and the number of filled pentagons, respectively. These values vary in a way that the total 
number of atoms in each fullerene group is preserved. Locations of X and Y on the fullerenes are chosen caring the 
maximum homogeneity or symmetry of their distribution. Using density functional theoretical (DFT) methods, we 
examine the relative stabilities of the fullerene groups by computing the relative energies, AE, the binding energies 
per atom, E^, NICS values and the band gaps, Eg. In doing this, we use a new hole density definition which is a 
generalization of the hole density for boron sheets 

hexagon hole number 

rjs — — 

Total number of atoms in original triangular sheet 

to spherical surfaces. This is done via Euler's law of polyhedra which states that a sphere can not be tessellated 
only by hexagons, but it also requires pentagons. Consequently, a naive definition of hole density for fullerenes would 
treat empty pentagons and as well as the empty hexagons as holes. However, this definition would be insensitive to 
capping procedure such that all members of the fullerene family ^^qXx+y would have the same hole density. For 
example, if the pentagons were treated as holes all the members of B^-^ family would have ?/=16/92. Our aim is to 
define the hole density so that each member in the family will have a different hole density depending on the values 
X and Y which affect the total energy. This, in turn, yields a non-trivial dependance between stability i.e energy 
and the geometric configuration. To accomplish this, we use the Euler's formula for polyhedron which restricts the 
number of pentagons to 12 and the number of hexagons to V/2 — 10. Here, V denotes the number of vertices or the 
number of atoms in our case. Then, the total number of atoms is given 'as V = 2(10 + H) where H is the number of 
hexagonal holes, and it is determined from the relation 60n^ ~ 2(H + 10). This yields the hole density definition for 
60n^ fullerenes as 



30n2 -IQ-X 

Henceforth we set n — 1, and drop the subscript. The lowest energy isomer of the each fullerene group and their 
geometric configurations are given in Fig. [T] Our analyses suggest a strong correlation between the relative energy 
and the cage geometry (See Fig. [2]). This is reflected in the fact that each lowest energy isomer in each fullerene group 
has the same normalized hole density ?7norm, which is defined as the difference between the fullerene hole density of the 
most stable cage, and the minimum hole density of the chosen fullerene family. In what follows, we define Ari[X, Y) 
which gives the difference between the maximum and the minimum hole densities in each fullerene family. Fig. [3] 
shows the distribution of the boron cage families with respect to Ary(X, Y) and the number fullerene members. It 
turns out that both distributions are gaussian and the fullerene families with the largest number of members has the 
largest Ari{X,Y). This fact will be of interest as the stability of the fullerene structures are concerned. In our case 
such fullerenes, constructed from Bgo skeleton, are given as Bg2, Bgo, Byg and Byg. 

The stability analysis requires the relative energies of boron clusters which depend on the choice of exchange- 
correlation functionals. Relative stability analysis of 2D and 3D boron isomers reveal that PBEO and TPSS functionals 
are more favorable over hybrid B3LYP [31 [51]. Specifically, Li et al. found that PBE, PBEO, TPSS and TPSSh 
basis sets are more reliable to describe the relative stabilities of Bgo isomers [3T]. Hence, all the DFT computations 
were performed using the Gaussian 09 [35j quantum chemistry package and Perdew-Burke-Ernzerhof (PBE) |36j and 
Tao-Perdew-Staroverov-Scuseria (TPSS) [37] functionals were used throughout the study along with the 6-31G(d; p) 
basis set. The harmonic vibrational frequencies for all the boron cages have also been calculated at the same levels 
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FIG. 2: (Color online) Relative energies (eV) as a function of fuUerene hole densities in the a) B76 b) Brg and c) B76 isomers 
calculated using PBE and TPSS density functionals and 6-31G(d; p) basis set and d) Bso isomers proposed and calculated by 
Li et al. [31] using five different density functionals and 6-31G* basis set. 



and all the stationary points were characterized as energy minima. NICS [351133] at the centers of the PBE/6-31G(d; 
p) optimized cages are computed at the gauge-independent atomic orbital (GIAO) |40H42] method with M06-2X/6- 
31G(d; p) level of the theory. 

Relative energies, symmetries, fuUerene hole densities, binding energies per atom, band gaps and NICS values of 
totaly 20 boron isomers are given in Table 1. X and Y values are included for all the studied boron cages. All the 
structures are fully optimized using PBE and TPSS density functionals with 6-31G(d; p) basis set. To compare the 
results, we plot the energies with respect to the hole densities for all the fuUerene groups shown in Fig. |2] Our results 
reveal that most stable isomers are B^g'^, B^f'^ and Blf'^. Energetically the lowest configuration of each fuUerene 
group shows a pattern such that hole density of the lowest energy isomer is shifted by 0.025 as the atomic number 
varies by 2 (see Fig. [2] and Table 1). When normalized, the hole density for the most stable fullerenes in each group 
approaches to the value rinorm= 0.07. We compare these results with the relative energies of Bgo isomers studied by 
Li et al. for five different density functionals shown in Fig. [2ji [3T]. Among them, the most stable isomer B|^^6 has the 
same normalized hole density, rinorm= 0-07, for all the basis sets. Numerical results suggest that the most stable cages 
are found in the neighboring groups of the Bgo group, which has the largest Ary(X, Y) value and the largest number 
of members (see Figure [s]). Theoretical simulations on less stable fullerenes like Bgo and B72 are in accordance with 
these results [M^- 

Electron delocalization or aromaticity plays an important role in energy lowering of molecules. NICS values can be 
used as a measure of aromaticity, where negative value means aromaticity, and positive value shows antiaromaticity. 
Computed NICS values, shown in Table 1, reveals that all the most stable fuUerene isomers B^g'^, B^l' ^ and Bgj'* are 
strongly aromatic with the values -66.4 ppm, -44.8 ppm and -24.7 ppm, respectively. To best of our knowledge, B^g'^ 
cage seems to have the largest negative NICS value -66.4 ppm with higher than any reported value for boron cages, 
including the core-shell structure Bioi, Bgo and B38 cages reported in [3TJ|43j. As an another stability measure, we 
compute HOMO-LUMO energy gap (Eg) [HHIS] and they are given in Table 1. The lowest energy isomers Byg-^, 
By|-^ and Bg|-^ have energy gaps 0.32 eV , 0.40 eV and 0.56 eV ^ respectively. The same trend for the relative energies 
of the lowest energy isomers also persists for the band gaps; the most stable cages have almost the largest Eg values. 
However, their values are low compared to the Bgp buckyball with appreciable energy gap of 0.99 eV . 

In Table 2, we give the binding energies per atom, the band gaps and NICS values of the well known a-BgQ-^ [15j . 
V-shaped Bg5^^ [17] and 7-Bfoo' [I^- In accordance with the previous results [3T], we found that Bgo buckyball has 
higher E;, than V-Bgo and they have antiaromatic NICS values which are given as 10.5 ppm and 22.9 ppm respectively. 
On the other hand, 7-B^q5° has a negative NICS value of -9.6 ppm but lower E^ which is in contrast to the result 
presented in jl2j . A comment related to the stability of Bioo family is in order. The 7-B100 might not be the most 
stable fuUerene structure of Bioo family and of all the remaining boron families which are constructed by capping the 
Bgo skeleton (adapted from Cgo). Because, as discussed in a preceding paragraph, the largest hole density separation, 
A?7(A', F), belongs to the Bhq isomers not Bioo- Therefore, the Bhq group and maybe other fuUerene families with 
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TABLE I: Symmetries, Relative energies, AE (eV), fuUerene hole densities, ri{X,Y), binding energies per atom, Eb (eV), 
NICS (ppm) and HOMO-LUMO energy gap. Eg (eV), values of boron isomers. NICS and Eg values are obtained at GIAO 
M06-2X/6-31G(d; p) and PBE/6-31G(d; p) levels, respectively. Relative and binding energy results in the first and second 
rows of each fuUerene isomer are calculated using FEE and TPSS density functionals, respectively. 
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FIG. 3: (Color online) a) Change of /S.r]{X,Y) with respect to fuUerene systems B6o-92. Largest separations are between 
B70-82. b) Change of the number of members of each fuUerene group with respect to fuUerene systems B6o-92. Largest number 
of members belong to B70-82 • 



larger A?7(X, y) than Bioo is expected contam energetically more stable member(s) with larger binding energy than 
that of 7-Bioo. This will be analyzed in more detail in a futm'e study. 



TABLE II: Symmetries, binding energies per atom Eb (eV), NICS (ppm) and HOMO-LUMO energy gap, Eg (eV) values. 
Binding energy results in the first and second rows of each fuUerene system are calculated using PBE and TPSS density 
functionals, respectively. 

Eb NICS Eg 
a-Big^°(Th) 5.77 10.5 0.99 
5.36 

V-Bio^2(Th) 5.76 22.9 0.32 
5.36 

7-B?!?5°(C2h) 5.75 -9.6 0.16 
5.35 
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In summary, we have carried out DFT calculations at the level of PBE and TPSS on 20 different boron isomers 
constructed originally from Bgo. Among them, the most stable isomers have the same normalized hole density, r]„orrm 
if the pentagons are not included in the definition of hole density eq{l] Distribution of boron cage families with respect 
to Ar]{X, Y), and number of members in each family is given. Our analysis shows that the most stable cages are found 
in the vicinity of the fuUerene group with the largest Ai]{X, Y) and the largest members along with the other fuUerene 
groups. In this respect, we predict the Buo family and other cages with large Arj{X, Y) as energetically more favorable 
over 7-3^05". Computed geometries do not cover all the possible fuUerene groups which preserve the total number of 
atoms, for example; Bs2 has also B|2'"^, B82'^, Bgj', Bg2'^ and Bg2'^^ isomers. Therefore, we think that the normalized 
hole density of Bgo+x+r molecules will be around 0.07. This study shows that boron fullerenes exhibit correlations in 
respect to stability and geometric structure as in the boron sheets. In this sense, we have introduced a new stability 
measure for boron cages which are constructed from the same core fuUerene; we evaluate the energies of the fullerenes 
and relate it to the hole density. These results might give insight into the evolution of boron stuructures; specifically 
in understanding the relative stabilities of B12 containing stuffed boron cages [511 ITftil^ and boron fullerenes. 
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